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INVESTIGATION  OF  IMPACT  RESPONSE  FOR  CFRP/STEEL  HYBRID  COMPOSITE 


The  response  of  composite  plate  embedded  with  steel  wires  subjected  to  low  velocity  impact  was  studied. 
The  steel  wires  were  embedded  within  the  layers  two  and  three  of  the  composite  laminate  plates.  The  first  order  shear 
deformation  theory  (FSDT)  as  well  as  the  Fourier  series  method  was  utilized  to  solve  the  governing  equations  analytically. 
A  computer  program  was  developed  based  on  the  analysis.  The  effect  of  volume  fraction  of  steel  on  the  deflections  and 
in-plane  strains  and  stresses  on  the  composite  plate  under  impact  was  studied.  It  was  seemed  that  the  embedding  the  steel 
wires  within  the  layers  of  laminated  composite  plates  would  improve  the  impact  resistance  of  the  plate  where,  the 
transverse  deflection  of  the  center  of  the  plate  reduced  from  (1.68E-07  m)  in  the  composite  medium 
(the  composite  without  steel  wires)  to  (1.25E-07  m)  in  the  hybrid  composite  plate  in  which  the  volume  fraction  of  the  steel 
is  9.3%. 
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INTRODUCTION 

Composite  materials  which  are  reinforced  with  fibers  have  many  advantages  due  to  their  high  specific  strength 
and  specific  stiffness.  This  can  be  fabricated  to  optimize  the  strength  and  stiffness  by  changing  fabrication  parameters  such 
as  the  stacking  sequence  and  the  number  of  plies  of  hybrid  materials. 

In  many  applications  such  as  aircraft,  automobiles  industries  and  sport  equipment  the  composite  plates  have  been 
used.  The  dynamic  response  of  composite  structures  subjected  to  dynamic  loading  has  been  studied  in  terms  of  analytical, 
numerical  [ ' 21  and  experimental  work  [ '  I  Theoretically,  many  researches  have  been  developed  for  studying  the  behavior 
of  composite  structures  under  low  velocity  impact.  The  analytical  model  presented  by  Pierson  and  Vaziri  [51  based  on  the 
combined  effects  of  shear  deformation,  rotary  inertia  and  the  nonlinear  Hertzian  contact  low  with  an  aim  of  studying  the 
low  velocity  impact  response  of  simply  supported  laminated  composite  plates. 

Embedding  the  Shape  Memory  Alloy  (SMA)  wires  inside  the  traditional  polymer  composite  can  reinforce  the 
damage  tolerance  of  multilayered  composite  structure.  This  is  because  the  existed  SMA  wires  can  generate  recovery  tensile 
stress  inside  the  structure  and  hence  reduce  the  deflection  and  the  in-plane  strains  and  stresses  in  transverse  loading  . 
The  effect  of  low  velocity  impact  on  the  composite  structures  and  prestressing  of  these  structures  was  extensively  studied 
in  the  past  by  Abrate  ,  Olsson  [8~91,  and  sun  et  al. [10]  and  many  others.  Birman  et  al.  also,  demonstrate  that  if  the  SMA 
wires  are  embedded  within  the  traditional  polymer  composite  plates,  tensile  stress  that  results  in  heating  of  the  wires 
increase  the  impact  resistance  of  the  structure. 
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The  constitutive  relationships  are  obtained  following  the  method  of  Pates  et  al. [12]  and  Zhong  et  al. [  '.  According 
to  this  method,  micro  mechanical  relationships  are  obtained  in  two  phases.  First,  these  equations  are 

Obtained  for  a  composite  medium  consisting  of  carbon  fiber  and  matrix.  Second,  these  relations  are  formulated 
for  a  composite  medium  treated  as  an  orthotropic  matrix  with  steel  embedded  within  a  medium. 

In  this  paper,  low  velocity  impact  of  symmetric  cross  ply  composite  plate  is  considered.  Steel  wires  are  embedded 
within  layers  in  xx  direction.  The  steel  effect  is  employed  to  reduce  transverse  deflection  of  the  plates  and  stresses. 

Assumptions 

Consider  a  cross-ply  composite  plate  subjected  to  low  velocity  impact  at  the  center.  Steel  wires  are  embedded 
within  the  layers  two  and  three  which  result  in  ahybrid  composite  plate.  The  analysis  is  based  on  a  number  of  assumptions 
that  are  listed  below: 

•  Perfect  bonding  exists  between  layers. 

•  Low-velocity  impact  can  be  represented  by  a  half  sine  function  of  time. 

•  The  carbon,  steel  and  matrix  experience  the  same  axial  strain  in  the  fiber  direction  and  the  same  transverse  stress 
in  the  plane  perpendicular  to  the  fiber. 

•  The  first  order  shear  deformation  theory  (FSDT)  as  well  as  the  Fourier  series  was  used  to  solve  the  governing 
equations  analytically. 

MICROMECHANICAL  MODELS 

The  method  employed  to  develop  a  micromechanics  of  steel  hybrid  composite  is  similar  to  that  used  by 
Pates  et  al.  [U]  and  Zhong  et  al.  [  I  According  to  this  method,  the  properties  of  a  composite  medium  that  consists  of 
ordinary  fiber  and  matrix  are  evaluated  first  using  an  appropriate  micromechanics  then  the  addition  of  steel  wires  is 
considered  using  the  same  micromechanical  model  with  composite  medium  treated  as  a  transversely  isotropic  matrix. 
The  formula  for  evaluation  of  the  properties  of  a  hybrid  steel  composite  ply  is  listed  in  the  APPENDIX. 

Constitutive  Equations 

The  plane  stress  constitutive  equations  for  a  lamina  of  constant  thickness  in  the  laminate  coordinates  are 


Where  Ql;the  transformed  reduced  stiffness  matrix,  {o~}and  {e}  are  the  stresses  and  strains  in  laminate  coordinates 
system  respectively.  Using  integration  of  equation  (1)  through  the  thickness  of  the  plate  (h)  yields 


(1) 


{QHksh  A 


]{y}ij=4,5 


(2) 
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Where,  N  and  Q  are  a  vector  of  forces  and  M  is  the  vector  of  moments.  Atj,  Btj,  and  Dtj  are  the  components  of 
extensional  and  shear,  coupling  and  bending  stiffness  matrices  respectively.  e°Andy  are  the  mid-plane  and  the  shear  strain 
respectively  and  k  is  the  curvature.  In  addition,  ksh  is  the  shear  correction  factor  and  is  often  taken  to  be  (n2/12). 

Governing  Equations 

The  assumed  displacement  field  is  [14l 

u(x,  y,  t)  =  u° {x,  y,  t)  +  zipx  {x,  y,  t) 

v  {x,  y,  t)  —  v°  (x,  y,  t)  +  zipy  [x,  y,  t) 

w  =  w°(x,y,t)  (3) 

Whereu0,  v°  and  w°are  the  plate  displacement  in  x,  y  and  z  direction  of  the  plate  mid-plane  and  xpx  and  ipy  are 
the  shear  rotation  in  the  x  and  y  direction. 

Reducing  the  equations  (2)  to  a  simply  supported  symmetric  cross  ply  laminated  plate, 
(By  =  A16-A26  —  D16  —  D26  —  A45  —  0).  Also  by  neglecting  the  rotary  inertia  (73),  results  the  following  equations. 


dy  dy2 


Dl1  -dx^+°12dx^  +  °66  Vdy^  +  dx^)'  KshAss  Vx  +^x~)-° 

d2ip 


dxdy 

°^  Kaxdy  +  dx*  )  +  °12  dxdy  +  ° 

Where 


(4) 


(Atj,  BtJ,  Dy)  =  S^2/2  Qij{l,z,z2)  dz 


(P.h.h)=S\Po  (XXz2)dz 


(5) 


Where  p  represent  the  density  of  the  plate,  p0  the  density  of  each  layer,  h  is  the  thickness  of  the  plate, 
Qijfor  i,j=l,  2,  6  are  the  reduced  in-plane  stiffness  matrix  components  and  Qij  for  i,j=4,5  are  the  reduced  transverse  shear 
stiffness  matrix  components. 

The  work  was  focused  upon  a  simply  supported  rectangular  plate  with  dimensions,  length  (a)  width(b)  and  with 
the  boundary  conditions: 


i//y=0atx=0,  a  also  w0  —  ipx  —  Oat  y=0,  b 


ANALYSIS  OF  THE  RESPONSE 

The  analysis  of  response  (stress,  strain  and  displacement  )  is  based  on  the  expansion  of  force,  displacement  and 
rotation  in  a  double  Fourier  series,  each  expression  is  assumed  to  separately  consisting  of  position  and  function  of  time, 
therefore,  the  above  boundary  condition  are  satisfied  the  following  expressions  (Navier  solution). 
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w0  (x,  y,  t)  =  £n=i  Sm=i  Wmn  00  sin  ax  cos  py  1  (6) 
^xO,y,t)  =  Sn=iSm=i^mn(t)cosaxsinPy 

xpy(x,y,t)  =  Xn=iSm=i^mn(t)sinaxcosPy 

Where,  the  Umn,  Vmn,  Wmn,  Xmn  and  Ymn  are  the  time  dependent  coefficient  to  be  determined  and  a  —  mn  /  a 
and  p  —  nn  /  b 

In  the  case  of  concentrated  force  F(t)  applied  at  the  location  (x,  y)  the  forcing  function  is  equal  to [15~161 


q(x,y,  t)  =  YimYinVmnit)  sm{mn  /  d)  x  sin  (nn/b) 


(7) 


pmn(t)  are  the  terms  of  the  Fourier  series.  It  is  clear  that  the  maximum  deflection  is  encountered  if  the  impact 
occurs  at  the  center  of  the  plate  where  (x  —  xc  —  a/2)  and  (  y  —  yc  —  b  /  2).  Therefore,  the  largest  deflection  and  stress 
will  be  occurred  also  at  the  same  location.  For  a  concentrated  force  located  at  the  point 


(*c,yc)imVmn(.t)  =  (AFCO/ab).  sin  (™)  xc.  sin  (™)yc 

Where  F(t)  is  the  impact  force  and  can  be  represented  [18] 
F(t)  =  p  sin(7rt/t0)ForO  <  t  <  t0 
F(t)  =  OFort  >  t0 

Where  to  is  the  impact  duration  and  pits  amplitude. 


(8) 


(9) 
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Using  equation  (6)  the  system  of  equations  (4)  can  be  reduced  to  the  following  system  of  ordinary  decouple 
differential  equations. 


(PmnOO) 


Where,  the  parameters  of  the  above  equation  are  defined  as  follow 

mn  —  m22  —  m33=i1>ku  —  kshiAss^xx  +  A^dyy),  k12  =  k]hA55dx ,  k13  =  kghA4idy 
k2i  =  ^12  '^22  =  kshA55  +  (fndM  +D66dyy)  ,  k23  —  (D12  +  D66)dxdy 

^31  =  ^13  '  ^32  =  ^23'  ^33=^sh^44+  iPe6^-xx   +  D22dyy) 

Based  on  the  analytical  model  and  using  matlab  software,  can  calculate  the  parameters  Wmn,  Xmn  and  Ymn  and  by 
substituting  these  values  into  (6),  the  values  ofw0  ,ipx,ipy  would  be  calculated.  Using  [  \ 


ex  =  Z-  $xx>ey  =  z-  ^y.y  •  Yxy  =       {$x,y  +  $y,x)  • 

NUMERICAL  CALCULATIONS  AND  DISCUSSIONS 


(11) 


By  using  the  computer  program  for  the  following  numerical  example  can  illustrate  the  effectiveness  of  the  steel 
wires  in  limiting  a  damage  inflicted  by  low-velocity  impact  where,  the  effect  of  embedding  the  steel  on  dynamic  response 
of  composite  plate  under  impact  is  studied  and  discussed.  Asymmetric  cross-ply  simply  supported  for  all  edges  rectangular 
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plates  made  of  carbon  fiber-polyester  with  different  volume  fraction  of  steel  and  with  a  lamination  of  [0°  90°  90°0°]  was 
considered.  The  size  of  the  plate  is  24  cm  X  24  cm  X  0.6  cm.  The  impact  characteristics  are  p  =  300  N,  to=300  jis  and 
volume  fraction  (  vs)  of  steel  was  taken  equal  to  (0%,  2.25%,  3.25%,  9.3%).  The  plate  materials  considered  in  this  work  as 
mentioned  were  carbon  fiber/  polyester  with  the  following  properties. 

•  Carbon  Fiber 

E1  =  230  GPa,E2  =  E3  =  15  GPa,G12  =  27  GPa,G23  =  G13  =  7  GPa,v12  =  0.2,  p  =  1790  Kg  /  m3 

•  Polyester 

E  =  3  GPa,  G  =  1.6  GPa,  v  =  0.35,  p  =  1100  Kg  /  m3 

And  the  steel  wires  of  (0.7  mm)  diameter  had  the  following  properties 

E  =  193  GPa,G  =  77  GPa.v  =  0.27, p  =  8000  Kg  / m3 

The  accuracy  of  present  analytical  and  used  program  was  verified  by  a  comparison  of  deflection -time  relationship 
obtained  from  this  solution  with  those  obtained  using  a  finite  element  solution  of  .  In  this  comparison  the  volume 
fraction  of  the  steel  was  taken  to  be  zero.  The  steel  wires  were  placed  along  the  x -direction,  which  corresponds  to  the  90° 
fiber  orientation  in  composite  medium  and  only  in  the  layers  2  and  3. 

Figure  1  shows  the  maximum  value  of  deflection  (w)  decreases  from  (1.68E-07)  in  the  composite  medium 
(the  composite  without  steel  wires,  curve  1)  to  (1.25E-07)  in  the  hybrid  steel  composite  plate  in  which  the  steel  have  the 
volume  fraction  of  9.3%,  curve  4.  Thus,  26%  reduction  is  occurred. 

x  10"7 
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-1 
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0  1  2  3          4  5  6 

Time  (sec)  x  1Q-4 

Figure  1:  Effect  of  Steel  Volume  Fraction  on  Deflection  History 

The  reduction  of  both  the  xx  and  yy  direction  of  the  in-plane  strains  and  stresses  during  the  impact  are 
remarkable  where,  maximum  stress  and  strain  reduced  by  8%,  4%,  28%  and  28%  for  stress  xx,  stress  yy,strain  xx  and 
strain  yy  respectively  in  the  composite  plates  to  hybrid  composite  plates.  Figures  2,  3,  4  and  5  demonstrate  the  variation  of 
the  in-plane  strains  and  stresses  with  time. 
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Figure  2:  Variation  of  Strain  xx  with  Time 
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Figure  4:  Variation  of  Stress  xx  with  Time 
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Figure  5:  Variation  of  Stress  yy  with  Time 


CONCLUSIONS 


The  impact  response  of  the  low  impact  velocity  upon  hybrid  steel  composite  plate  is  studied  using  simultaneously 
the  first  order  shear  deformation  theory  (FSDT)  and  Fourier  series  to  solve  the  governing  equations  analytically. 
A  computer  program  based  on  the  analytical  model  and  using  the  matlab  software  was  utilized.  The  program  is  capable  to 
compote  the  following  parameters: 

•  Deflection  history  of  the  composite  plate 

•  In-plane  strains  and  stresses  history  of  the  composite  plate 

The  results  obtained  from  this  program  also  indicated  that  the  embedding  the  steel  wires  within  the  layers  of  carbon 
fiber-epoxy  laminated  composite  plate  results  in  that  the  deflection  and  in-plane  strains  and  stresses  of  the  plate  decrease 
with  the  increasing  the  volume  fraction  of  steel.  Furthermore,  the  effectiveness  of  steel  wires  can  be  further  improved  by 
optimizing  their  distribution  through-out  the  plate. 
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APPENDICES 


Micromechanics  of  steel  hybrid  ply 


Properties  of  a  composite  medium  consisting  of  transver  selyiso  tropic  carbon  fiber  and  an  isometric  matrix 
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Properties  of  a  hybrid  steel  ply 


El  =  E{vc  +Esvs  ,  El  = 


n-v^-t)] 


.  V12  =  vf2Vc  +  VSVS 


In  these  equations,  the  superscript,  'c'  identifies  the  composite  medium  (fiber  and  matrix),  and  'h'  identifies  the 
hybrid  medium  (composite  medium  and  steel),  while  the  subscripts  'm'  and  '/'  refer  to  the  matrix  and  the  fibers, 
respectively.  The  subscript  T  and  '2'  refer  to  the  longitudinal  along  the  (fibers)  and  transverse  directions  respectively. 
Other  notations  are  as  follow: 

E=  modulus  of  elasticity,  G=shear  modulus,  v=ppoison's  ratio,  v=  volume  fraction,  where  the  subscript  Y 
identifies  the  steel. 
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